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Inclusive S- and P-wave charmonium productions in the bottomonium ground state r]i, decay are 
calculated at the leading order in the strong coupling constant ctg and quarkonium internal relative 
velocity v in the framework of the NRQCD factorization approach. We find the contribution of 
'?!) ~^ Xcj + 5 5 followed by Xcj ~* J/ti' + 'y is also very important to the inclusive J/^p production in 
the rjb decays, which maybe helpful to the investigation of the color-octet mechanism in the inclusive 
J/i/' production in the r/j, decays in the forthcoming LHCb and SuperB. As a complementary work, 
we also study the inclusive production of ric, and Xcj in the rib decays, which may help us understand 
the X(3940) and X(3872) states. 



PACS numbers: 12.38.Bx, 12.39.Jh, 13.20.Gd 



I. INTRODUCTION 



The existence of the spin-singlet state rj^, which is the 
ground state of hh system, is a solid prediction of the 
non-relativistic quark model. Since the discovery of its 
spin-triblet partner T, people have make great efforts to 
search for it in various experimental environments, such 
as in e+e~ collisions at CLEO in 77 collisions at LEP 
II Hand in pp collisions at TevatronQ. Unfortunately, 
no evident signal was seen in these attempts. Recently, 
a significant progress has been achieved by Babar collab- 
oration. After analysing about 10^ data, they observed 
77b in the photon spectrum of T(3S') — > 7?7b[3| with a 
signal of 10 ij significance. They found the hyperfine 
T(IS') - i]b mass splitting is 71.4j:|5(stat) ± 2.7(sjyst) 
MeV. Soon after, it was also seen in T{2S) JilbM by 
another group in Babar, and the mass splitting is deter- 
mined to be 67.4t:^;^(stat) ± 2.0(syst) MeV. 

On the theoretical side, considerable works have been 
done to study its properties. The mass of rji, has been 
predicted by potential model effective theory]^ and 
Lattice QCD \§\. And the recent determinations of 
T(LS) —Jib mass splitting in the range of 40 — 60 MeV 
0, Eol EH, are consistent with the Babar's results. 
Aside from its mass, the production and decay proper- 
ties of r]i, have also been considered. The number of rjh 
produced in e~^e~ ^ j + rjb at S-factoriesfTs] is found 
to exceed that produced at LEP II by about an order of 
magnitude. In Ref.[l3], the authors calculated the pro- 
duction rates of rji, at Tevatron Run II and suggested 
to detect it through the decay of rji, — > J/4>J/4'i while 
in Ref . . it was thought that the double J/ip chan- 
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nel might be overestimated and it was suggested the 
rjb D*D^*^ channel to be the most promising chan- 
nels. An explicit calculation of rjb — > J/tpJ/tp at NLO 
in ti^[l^ and NLO in as[l3 shown that this branching 
fraction is of 10^^ order which is about four orders of 
magnitude smaller than that given in Ref.[i|. And the 
author in Ref. argued the effect of final state interac- 
tions in rib DD* —^ J/ipJ/ip was also important. Some 
other exclusive decay modes such as rib 7«//'00) US] 
and rib decays into double charmoniaj2l|, and inclusive 
decays, e.g. rib ccccfTsj and rib J/i^ + X[22] have 
also been taken into account. 

However, comparing to the cc ^Sq state 77c, our knowl- 
edge about rib is quite limited, and doing some further 
works is necessary. In this paper, we will systemically 
study the inclusive decays of rib into S- and P- wave char- 
monium states. The motivations of this work arc four- 
fold. First, in these processes, the typical energy scale rrib 
in the initial state and TOc in the final state are both much 
larger than the QCD scale Aqcd , so we can calculate the 
decay widths perturbatively and the non-perturbative ef- 
fect plays a minor role, which will reduce the theoretical 
uncertainties. Second, the branching fraction of the in- 
clusive decay process is much larger than that of the ex- 
clusive process, which makes the test of theoretical pre- 
diction for the inclusive process be more feasible. Third, 
in Ref. [22], Hao et al. have calculated the branching ratio 
of rib — > J /tp + V and found the contribution of the color 
octet process rib cc('^S'[*') + g is larger than the one 
of the color singlet process by about an order. Since the 
color-octet process also contributes to P-wave states XcJ 
production, in which the Xci and Xc2 has about 36% and 
20% branching ratio to J/ip-\-^ respectively, so we expect 
that the contribution of rib ~* Xc.i + V process followed 
by XcJ ~* J 1^ + 7 might also be important for inclu- 
sive J /ip production in rib decay. Fourth, in recent years, 
many charmonium or charmonium-like states have been 
found at i3-factory (see Ref. p3l [23. [25| for a review). In 
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the further coming LHCb and Super-B, when accumu- 
lating enough data, it might be possible to observe the 
interesting decays of r]b to X(3940) or X(3872) etc.. 

The J/ip inclusive production has already been 
studied[24|, and the J/ip{r]c,Xcj) production in associ- 
ation with cc pair has been discussed in our previous 
workpsj. As important supplements, here we are going 
to consider the contribution of rji, ??c(Xcj) + 99 process 
in the non-relativistic limit at leading order in 0:3. 



II. NRQCD FACTORIZATION FORMULISM 

Because of the non-relativistic nature of bb and cc sys- 
tems, we ado pt t he non-relativistic QCD (NRQCD) ef- 
fective theory [231 to calculate the inclusive decay widths 
of rjb to charmonium states. In NRQCD, the inclusive 



decay and production of heavy quarkonium are factor- 
ized into the production of short distance coefficient and 
the corresponding long distance distance matrix element. 
The short distance coefficient can be calculated perturba- 
tively through the expansion of the QCD coupling con- 
stant as- The non-perturbative matrix element, which 
describes the possibility of the QQ pair transforming into 
the bound state, is weighted by the relative velocity vq 
of the heavy quarks in the heavy meson rest frame. 

In the framework of NRQCD, at leading order in vi, 
and Vc, for the S'-wave heavy quarkonium production and 
decay, only the the QQ pair in color-singlet contributes. 
For P-wave Xcj production, the color singlet P-wave ma- 
trix elements and color-octet S'-wave matrix element are 
both in the same order of v^. Then the factorization for- 
mulas for the processes under consideration in this work 
are given by: 



(la) 



XcJ + X) = fi(65(isW) ^ cc(-VW) +X)(ry,|a(isW)h,)(0?-^(3p[il)) 



(lb) 



short-distance factors 



and 
and 



where the Fs are the 

{O^"-' sf'^)) are the long-distance matrix elements. 
During our calculation of the short distance coefficients 
associating with the P-wave color-singlet matrix el- 
ements, there will appear infrared divergence. This 
divergence will be absorbed into the color octet matrix 
element {O^"-' {^sf)). 
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FIG. 1: Typical Feynman diagrams for the short distance pro- 
cess: (a) bbl'-So, 1] cc[^4''(^-P™)]+ff5; and (b) bfef^So. 1] 



III. rib^ric + gg 

We first consider the S-wave rjc production from rji, 
decay. At leading order in as, there are eight Feynman 
diagrams for (66(^4^') ^ ccC-S^^^) + 99. The typical 
one is shown in Fig. [la]. The general form of the short 
distance coefficient can be expressed as: 

f((66(i5W) - cc-(isW) +gg) = ^f{r), (2) 

where r = mdmi, is a dimensionless parameter. Since 
there is no infrared divergence, we calculate /(r) dir ectly 
using the standard covariant projection technique [28|. 
Given mb ~ 4.65GeV, TOc = l.SGev, we get /(r) = 23.1. 
In NRQCD, up to order, the relations between the 
color singlet matrix elements and the non-relativistic 



wave functions are 

(0c^^C4")> = ^\Rism\^+o{vt)). (3) 

In order to compare with our previous work before, we 
choose the same numerical values with m,\, = 4.65GeV, 
nic = 1.5GeV,as = 0.22, |PJs(0)|2 = 0.81GeV^ and 
|p6^(0)|2 = 6.477GeV3[2^. Then we get 

r(77fc ^ Vc + 99) = 0.83 kev. (4) 



^ For the color-singlet four-fermion operators, there is a additional 
factor compared to those in Ref. |27|| 
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FIG. 2: The scaled energy distribution of rjc for the rji, 
rjc + gg process. 



The total width of rjb is estimated by using the two 
gluon decay, which at leading order in as and Vb is read 
to be: 



TTotai ^T{r]i,^ gg) 



2ai 



3m 



2 l^isl 



9.67MeV. (5) 



In our previous work, we got T{rii, ^ tJc + cc) ~ 
0.27keV[2y]. So the branching ratio of inclusive decay 
of ?7f, into is 



Br(77h -^ric + X) = 1.1 x 10" 



(6) 



in which the contribution of gg process is about 3 times 
larger than that of the cc process. The re-scaled energy 
distribution curve dT/dxi for rji, r]^ + X is shown in 
Fig. [2], where xi is the ratio of rjc energy E"^^ to mt. 

Recently the X(3940) state was observed by the Belle 
Collaboration in the recoiling spectrum of J/V' in e~^e~ 
annihilation 30] . It is most likely to be a ric{3S) statefsij. 
In the non-relativistic limit, the only difference between 
77c and ric{3S) is the value of wave function. If X(3940) 
is the ric{3S) state, we predict the branching ratio of 
X(3940) production in rj^ decay to be 



Br(77h ^ X(3940) + X) ~ 0.62 x 10" 



(7) 



To obtain the prediction, we have chose |-R35(0)P = 
0.455GeV^[2^ to take the place of |i?Js(0)|2 = O.SlGeV^ 



IV. 



Xcj + gg 



of which is shown in Fig. [lb]. Using the D dimension 
spin projector expression[32], at leading order in a^, the 
short distance factor is given by 



t{bbCSo) 
(47ra 



24m^r' 



^ $2 — 



c-c{^ Sf^)+g) = 
2}{D-3) 



(D-l) 



(8) 



r(l-e)(l-r^) 
87rr(2-2e) 



is the 2-body phase space 



where $2 = i:^} 
in Z?-dimensions. 

The calculation of the color-singlet coefficient in full 
QCD is a little more complicate. The Feynman diagrams 

for 66(i4^1)(P) ^ cc(3pW)(pi) + 3(^2)5(^3) is the same 
as those for rjc production process. Such 1^3 process 
can be described by the following invariants: 



2P^ 
M2 



(9) 



where AI = 2mi,. In D = 4 — 2e dimensions, the three- 
body phase space is given by 

d$(3) = KK{{aa + bb ~ Xi)(x2 + aa - bb))~^ 

(1 + — x\)^'^8(2 ~ xi ~ X2 — x^)dxidx2dx^llQ) 



where r — mc/mi,,aa = xf ~ Ar'^ /2, bb = (2 — Xi)/2 

and KK = 32^3r(2-2e) ■ 

In calculating the amplitude, we put the four dia- 
grams with the gluon carrying a momentum of p2 emit- 
ted from the charm quark line together and label their 
total amplitude with M2 and the total amplitude of 
the left four diagrams with the gluon carrying a mo- 
mentum of p3 emitted from charm quark line is repre- 
sented by M3. The total amplitude M = M2 + M3 and 
|M|2 = IM2P + IMsp + 2Re(MlM3). 

As being illustrated in Ref . [2 7| . for P-wave case when 
Pi{i — 2,3) goes to zero, there will be singularities in 
Ali. However, because of the four- momentum conserva- 
tion, P2 and P3 can not be soft simultaneously in the 
phase space. Therefore, the integration of the interfer- 
ence term 2Re(M2M3) is finite. We could perform it in 
4-dimensions directly. For the symmetry of the two glu- 
ons, the result of phase space integration for IM2P and 
|Af3p are equal to each other. We only need to calculate 
one of them. The total Fi then could be written as 



Fi = 2F 



Ah 



F 



Int- 



(11) 



As mentioned above, the color singlet short distance 
coefficients are infrared divergent in full QCD calculation. 
We will adopt the dimensional regularization scheme to 
regularize the divergence. To absorb the divergence into 
the color-octet matrix elements {O^"-' (^sf^)), we are {p2~P 
necessary to calculate the color-octet short distance co- 
efficient in Z? = 4 — 2e dimensions. The 66(^S'q^' 

-|- g process includes two Feynman diagram, one 



where Fm2 and Fjnt are the contribution related to IM2P 
and 2Re(M2 AZ^s) respectively. 

We now present how we calculate Fjv/a in detail. The 
denominator of charm-quark propagator in Fig. [la] is 



-2p2-Pi 



9c =0 



oc 



(l + r^-xi-xa), (12) 



where pc = ^ — gc is the momentum of anti-charm 
quark and qc is the relative momentum of c and c. 
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When cc in P-wave configuration, we need to know 
the first derivative of the amplitude with respect to Qc- 
Then in the non-relativistic Umit, there will be three 
kinds of the divergences in IM2P, which are proportional 

1 

to (i+r2_Ji_a;2)4, (i+r'i-xi-x2)i (i+r'i-xi-x2) ■ These 
terms will be divergent at point {xi,X2) = (1 + r^,0) 
which are not easily to be integrated out. We introduce 
two new variables (x'i,x'2), defined by 



Where f is finite and can be calculate in £) = 4 di- 
mensions. The phase space integration of the first term 
in Eq.lfTSl) is expressed as 



W (1 +r2 - xi)^ 



where 



l+r^ 



2r 



(1 



(17) 



X — X ^ Xi^ — 1 



1 + — xi 

X2 



(13) 



In the variables x'l and x'2 , the phase space is re-expressed 



as: 



l-{bb+aa) (1 ~ 



((aa + bb — X2){ 



1 — x'2 bb + aa 



))-'dx'2. (18) 



^2e^2-4. .l+r- .l-ibb-aa) 



327r3r(2 - 2e) 



2r 



-{bb+aa) (1 " ^2)'^ 



Furthermore, the integrals in Eq. (ll7p can be written into 
the sum of two terms defined by: 



(1 + -xi)'-^'{{aa + bb -X2)i 



1 — x'-, bb + aa 



(14) 



where aa = ^^^^ '^^ ' = ^2 — ^^jir^p^- 

And the three divergence structures are changed to be the 

fnrm nf 1 (1-^2)^ 1 (1-^2)^ „„A 1 (1-^2)^ 



respectively, which are all proportional to 
Then |M2p could be expanded as 



1 

(l+r2_j,.j2 



dx^ig(a;^i,6) ^ r'+^- dx',g{l + r\e) 
2r (l-Hr2-a;i)i+2^ (l-Hr2-a;;)i+2^ 



da:;(5(x;,e)-.g(l + r2,e)) 



(19) 



J2r (H-r2-x'i)l+2^ 

The first term on the right side includes - pole, and the 
second term is finite. Therefore we only need to keep the 
0{e) contribution when calculating g{l -t- r2,e) and the 
second term can evaluated directly by setting e = 0. 
Putting Eq.(ll) and (16) into together, we get 



/i(l + r2,4,e) 
(1 r2 - xiY 



+ /2(x;,4,e). (15) 



Accordingly, 



7^ T^div j f^fin 

^ M2 — ^ Mn + Mn 



(16) 



ri = 2(ni^ + rS?)-i-ri„t 



(20) 



calculated analytically, and f and f int are cal- 
culated numerically. For J ~ 0, 1,2, the expressions for 



r^}- are 



(J = 0) (21a) 



81 mlr^e 

MCACpT^^atKK (4 - Ar^ + 24 (l - r"^ (3 - 3r2 + r")) log(l - r^) + 12 (l + r^) log(r)) 

243m-^r5(-l + r2)^ 

^'^ ~ 81 m9 r5 e + 243m9r5(-l + r^f 

(2 - + 9r^ - 2r^ -f 3 (2 - Zr^ - 3r^ + 2r^) log(r) -I- 12 (l - Zr"^ + 3r'* - r^) log(l - r^)) [J = 1) (21b) 

fdiv _ 128 (-1 + r2) CACF{as^T^i'''YKK , U^CACpT^^aiKK 
1 i\,f„ — 



81m9r5e 1215m9r5 (-1 + r2)2 

(10 - 27r2 + 27r'' - lOr^ + 3 (lO - - gr"* + lOr^) log(r) - 60(-l + r^)^ log(l - r^)) (J = 2). (21c) 



The Ca — 3 and Cf = 4/3 in above equations are the color factors. It can be seen that for different J the diver- 
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gence part of f are the same, which wih be absorbed 
into the color-octet matrix element. And 21"^"^ +rint are 

^,4 



Olefin 

Mo 



ri„ 



.Aj{r) (forJ = 0,l,2). 



(22) 



When r = 1.5/4.65, we obtain Aa{r) ~ -9.71 x 10^, 
Ai{r) ~ -2.66 X 10^ and A2{r) ~ -6.06 x 10^. 

To cancel the infrared divergence of F^J^ , we also need 

to take into account the renormalization of {O^"-' ("^S"!^')). 
In MS scheme, it is given byjl^, fs^ 

^(i + log 4. - lE){^rj:{OreP'}')). (23) 



Combining the results in Eq. (8, 21, 22, 23), we finally ob- 
tain the infrared safe expressions for inclusive decay of rji, 
into Xcj{J — 0, 1, 2) states 



where Vi is 



27r2a3(l 



and rf are 



(24) 



'-{ri,\0,Cs'i%,){0^,-'{^sf^)), (25) 



^ ^-^<-^^'^;ff^)';^^)<^?:-(^^o')) (12(.'^ + 1) log(.) + 24(1 - 3.^ + 3.^ - r^) log(l - r^) + 
243TO^r''(l — r'^Y 

2(1 - r2)((61og2 - 5)r4 - 4(3 log 2 - -f 61og2 - 5 + 6(1 - r'^f \og{ — )) + 243^^(1 - ^')'-4o(r) ^^ea) 

MA Svr 



rl = ^"^".W^l^^l r ° ^^ (3(2^^ - 3r4 - ir^ + 2) log(r) + 12(1 - log(l - r^) + 



243m^r5(l - r'^ f 

+ (1 - r2) ((6 log 2 - 5)r4 + (7 - 12 log 2y + 6 log 2 - 5 + 6(1 - r^)^ log(^)) + 243r5(l - r^)2Ai(r) ^^eb) 

/iA 167r 



1215m^r5(l -r2)2 



(3(10r'' - 9r^ - Or^ 10) logr -f 60(1 - H)^ log(l - r^) + 



(1 - r^)(5(61og2 - 5)r* + (53 - 601og2)r^ + 5(61og2 - 5) + 30(1 - r^)'^\og{ — )) 

MA 



1215r5(l - r^^Aiir) 
Ifor 



)• 

(26c) 



It can be seen that the contribution of P-wave color- 
singlet is dependent on the factorization scale ma- When 
combining it with the color-octet 5— wave contribution, 
in which the matrix element also depends on /iA, the /iA- 
dependence will be canceled. 

To give numerical predictions, we also need to know the 
values of the long-distance matrix elements. The color 
octet matrix elements can be studied in lattice simula- 
tions, fitted to experimental data phenomenologically or 
determined through some other non-perturbative ways. 
Here we determined their numerical values with the help 
of operator evolution equations. In the decay process, 
the solution of the operator evolution equations are [13]: 



{Xcj\OsCSi;fiA)\Xcj) = (Xc,7|08('5i;/iAo)IXc,7> 

■In — -{Xcj\Oi{-'Pj)\xcj), 



(27) 



where (3q — ^^^^q "^^^ ■ We then naively relate the matrix 
element of production operator and that of the decay 
operator 0„ using 



{0^)^2J+l{H\0n\H) 



(28) 



When /iA 3> ^J■Ao^ the evolution term will be dominated, 
and the contribution of the initial matrix elements can be 
neglected. Since the operator evolution hold only down 
to scale rricV, we set the lower bound /iAo = and 
choose — 0.3. And we set /iA = 2toc since the di- 
vergence comes from the soft gluons linked with the cc 
pair. The P-wave color singlet matrix elements can be 
estimated through their relations with the first derivative 
of the non-relativistic wave function at origin which, in 
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non-relativistic limit, is given by 

(O?-^(M'')>«^^^^|i?^(0)p. (29) 

Setting Nf = 3, Aqcd = 390MeV and \R'M\^ = 
0.075 GeV^lii], we obtain 

r(?76 -> XcJ + gg) = (0.17, 1.55, 1.76)keV (for J = 0. 1, 2). 

The rjb XcJ + cc processes have been considered in our 
previous work, in which both the color-singlet and color- 
octet contributions have been included but with different 
values of the color-octet matrix elements [1^. If we use 
the color-octet matrix elements determined in this work, 
the results now become: 

^{Vb XcJ + cc) = 

(4.54, 4.21, 4.28) x lO^^^cV (forJ = 0, 1, 2), (31) 

which are about an order of magnitude less than the 
widths of rjb XcJ + gg processes respectively. Includ- 
ing the contribution of the associate processes, we then 
predict that the branching ratios for rjt inclusive decay 
into XcJ are 

Br(77fc ^ XcJ + X) = 

(0.22, 1.65, 1.87) X 10^^ (for J = 0, 1, 2). (32) 

The X(3872) state was discovered in pp collisions at 
Tevatron^SS'l and B decay at Belle [33|. Until now, people 
have not found an convincing explanation about it yet. 
The authors in [sH] suggest it is a Xci(2P) state. If it is 
a Xci(2P) state, we roughly predict 

Br(7?6 ^ X(3872) + X) = 2.25 x 10~^ (33) 

where we have chose |i?c(0)P = 0.102GeV^ and assumed 
the ratio between color-singlet and color-octet matrix el- 
ements does not change for 2P state. 

jo^-i'sf^)) ^ (gf^""^(^^f')) ,,,, 

(0?-(3pfl)) (Of(3«^^)(3Ffl))- ^ ' 

In [l^l, the authors have studied the ry?, J/ip+X pro- 
cess with r(77,, J/ip -f X) = 2.29keV. They found the 
contribution of color-octet process r/t, — > t//V'coior-octet + 
X is more than one order of magnitude larger than that 
of the color-singlet contribution. Since Xci and Xc2 could 



also decay to J/ip + 7 with Br(xci J/"^ + 1) = 36% 
and Br(xc2 J/'fp + 1) = 20%. The branching ratio 
of Xco J/i^ + 7 is so small that the contribution of 
this process can be neglected. Then re-scaling our result 
by the values of parameters in Reffl^l, we find the XcJ 
feed-down contribution to the decay of rji, into J/tp is: 

r(r/b ^ [J/t/j + 7)^^ J +X)^ 0.71keV, (35) 

which is about three times larger than that of color- 
singlet process. Therefore in the future experiment, when 
measuring the J/^ production in r]{, decay, the contribu- 
tion of rjb decays into XcJ followed by Xcj — *■ J/V' + 7 are 
also important. 

V. SUMMARY 

In this work, we have studied the inclusive produc- 
tion of charmonium state rjcXcj in the decay of ground 
bottomnium state rjb within the framework of NRQCD 
factorization formula. We find for the P-wave states Xcj 
case, the color-singlet processes 65(^5q^') cc(^P^^)+gg 
include infrared divergence. We show that such diver- 
gence can be absorbed into the 5- wave color-octet matrix 
element. To give numerical predictions, we use the poten- 
tial model results to determine the color-singlet matrix 
elements and estimate the color-singlet matrix elements 
with the help of operator evolution equations naively. We 
find that the branching ratios of r]b decays into r]c or Xcj 
plus anything are all of 10""* order. Furthermore, we 
also give the branching ratios of rjb X(3940) -I- X and 
rjb X(3872) + X, if the X(3940) and X(3872) are 
the excited 77c(35) and Xci(2P) states respectively. In 
Ref. ;22j. the authors investigated the color-octet mech- 
anism for J/^p production in rjb decay, our results show 
that the J/ip production from Xcj feed-down is also im- 
portant, since it is about three times larger than the di- 
rect J/ip production via color-singlet channel. These the- 
oretical predictions may not be observed in experiment 
for the time being, but will be very helpful to study the 
rjbS properties in the future experiment such as Super-B. 
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